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Abstract. We study the spherical model of a ferromagnet in d-dimensional cubes 
fin of volume |f2„| = n'^ and investigate large deviations of the magnetization of 
various domains C fi„. We focus our attention on the low-temperature regime, 
T < Tc, and consider domains Dk of three types: {d — l)-dimensional layers of 
width k, {d — 2)-dimensional rods, and KadanofF blocks. In the case of layers 
the large-deviation probabilities decay exponentially with n'*^^, and we obtain 
an explicit expression for the corresponding rate function. When the layer width 
fc <C the large-deviation probabilities are virtually independent of k. In the case 
of rods the probabilities of large deviations exhibit similar exponential decay, but 
this time it is distorted by logn corrections. In the case of Kadanoff blocks of size 
k the large-deviation probabilities decay exponentially with fc''^^. 

KEY WORDS: Critical phenomena; equivalence classes; correlation length; Kadanoff 
blocks. 



1 Introduction. 

Most models studied within the theory of critical phenomena describe the behaviour 
of various order parameters reasonably well. Qualitatively predictions of mean- 
field and of short-range finite-dimensional models are very similar in this respect. 
The only significant discrepancy is in the values of critical exponents. The latter 
predicted by 2D and/or 3D models differ quite substantially from the set of mean- 
field critical exponents. 

To find qualitative differences in the behaviour of various classes of models in 
the entire low-temperature region (below the critical point) one can look at the 
probabilities of large deviations: 

Pr [rriN e [a, b]] , 

where m^r is the corresponding order parameter, N is the number of microscopic 
degrees of freedom in the system, and (if we are indeed talking about large devi- 
ations) the interval [a,b] does not contain the equilibrium value, m*, of the order 



^Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, Dubna 141980, Russia 
e-mail: patrick@theor.jinr.ru 



1 



parameter. Typically, as N ^ oo, large-deviation probabilities exhibit the following 
asymptotic behaviour: 

— — lnPr[m7v € [a,b]] ~ min R{x), 

N x<^la,b] 

where R{x) > is the corresponding rate function. 

Although, being the probabilities of (extremely) rare events, the above quanti- 
ties may seem unimportant they have significant implications on, for instance, the 
existence/absence of such a mysterious phenomenon as the hysteresis loop. Large- 
deviation probabilities also describe the behaviour of thermodynamic systems when 
it is necessary to take into account some conservation laws. A frequently encoun- 
tered example is the conservation of the number of particles. Properties of typical 
configurations realising appropriate large deviations in a system where the number 
of particles is not conserved are closely related to equilibrium properties of the same 
system with the conservation law. 

For short-range finite-dimensional models the rate functions R{x) are always 
convex functions with a continuous first derivative. On the contrary, in the low- 
temperature regions the rate functions of mean-field models are not convex. There- 
fore, in this respect the properties of mean-field models are unphysical. 

For instance, the Hamiltonian of the Curie- Weiss model of a ferromagnet — the 
canonical mean-field model — is given by 




TT-CW 

-Hat = 

where Sj = ±1, for j = 1, 2, . . . , A^. A straightforward calculation, see, e.g., 
yields the following asymptotics 



— In Pr 

N 



If =A 




for k G {— A^, — A^ -|- 2, . . . , A^ — 2, A^}, where Pr[-] is the Gibbs distribution corre- 
sponding to the Hamiltonian if™ and the inverse temperature (3, 

F{x) = _ ^hx + ln(l -x) + ln(l + x), 

and k* minimizes F{-^). li h = and f3J > 1 the corresponding rate function 

Rc„{x) = F{x) - min F{y) (1) 

is not convex, see Fig. 1. When h the function Rcvj{x) has exactly one global 
minimum. If h (7^ 0) is sufficiently small i?cw(a^) also has a local minimum, which 
is often interpreted as a quasi- stationary state giving rise to the hysteresis phe- 
nomenon. 



2 





0.2 ■ 


. Rcw{m) 






0.1 ■ 






-1 ^ 






1 




m 




Figure 1: The rate functions for the magnetization of the Curie- Weiss model at 
13 = Pc = l/J (left) and at /3 = 1.2/ J (right). 



It is possible to consider a version of the Curie- Weiss model with continuous 
random variables Sj. In this case the corresponding rate function will differ from 
Rcwix) given by Eq. ([1]), but certain qualitative features of R{x) (for instance, the 
non-convex shape in the low-temperature region) will remain unchanged. Thus, al- 
though the rate functions R{x) are not universal quantities, some of their properties 
are identical within large classes of models. 

Short-range finite-dimensional lattice models can be divided into two broad 
classes: continuous and discrete. A well-known exactly solvable representative of 
the former class is the spherical model of a ferromagnet defined in cubical domains 
Qn of a square lattice Z'^, see [1] and the definition in Section 2. A cube fin con- 
tains n'^ = N lattice sites. The order parameter here is the magnetization rriN 
and the rate function describing the large-deviation probabilities for ttitv is given by 
Rsphix) = G{x) — G(0), where G{x) = ma.Xz>dg{x, z) and 



g{x, z) = — /3J(1 — x'^){z 



^^+2(2^/0 •••yo^^^---^^^^H' 



cos 



see [6]. The function R^^^{x) is always convex, see Fig. 2. In the low-temperature 
region P > (3^ this function vanishes on the entire interval x G [— m*,m*], where 



m* = y 1 — l3c/ 13 is the spontaneous magnetization obtained by switching off a ho- 
mogeneous external magnetic field h ^0. Switching on an arbitrarily weak magnetic 
field immediately leads to a strictly convex rate function with a positive second 
derivative and a unique minimum. Thus, when h ^ the spherical model on a 
finite-dimensional lattice does not have quasi- stationary states, and hence it does 
not exhibit any hysteresis phenomena associated with such states. 

Since the rate function Rsph{x) is equal to zero on the entire interval [— m*,m*] 
it is still necessary to find the first non-vanishing term in the large-n asymptotic 
expansion of InPr [rriN E [a, b]] when [a, b] fl [— m*, m*] 7^ . It turns out that in the 
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Figure 2: The rate functions for the magnetization of the spherical model at /3 = /3c 
(left) and at = (right). 



case of the spherical model with periodic boundary conditions this term is given by 



InPr [niN € [a, 6]] ~ —ri 



min i?S(x), 



where 



-^sph(-^) 



, for X E [—m* 



m 



see Fig. 3. Note that, unlike R{x) — the rate function of the order n*^, lower-order 
rate functions (the order n*^"^ here and the order n'^~^ below) do not have to be 
convex functions with a continuous first derivative. 

Among discrete models the 2D Ising model on a square n x n lattice is the most 
frequently studied example. The order parameter here is the magnetization rriN, 
N = 71^ . An explicit expression for the corresponding rate function, -Rising ('^), is 
not known. However, it is known that -Rising ('^) is always convex and, thus, it has 
the same shape as the analogous rate function within the spherical model shown 
in Fig. 2. In particular, for P > the function -Rising ("^) vanishes on the entire 
interval [— m*,m*], where m* = [1 — sinh~^(2/3J)]^''^ is the spontaneous magnetiza- 
tion. If [a,b] n [— m*,m*] ^ 0, then the leading asymptotics of the large-deviation 
probabilities is given by 



InPr [niiy E [a, b]] 



-n mm R^ (x), 

xG[a,b\ 



see [3l[?j|. It is quite remarkable that despite -Rising (a:^) is still unknown, in many cases 



one can find an explicit expression for the rate function r[1^j^^{x). For instance, it 
was shown in the paper |fB] that in the case of the 2D Ising model with periodic 
boundary conditions this rate function (in the notations adopted to ours) is given 
by 




for mo < \x\ < m*, 
for |x| < mo. 
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Figure 3: The rate functions i?gph(m) and -Rising("^) for ths magnetization of the 3D 
spherical and 2D Ising models in the low-temperature region (3 > (3c- 

see Fig. 3, where w is the surface tension associated with the droplet boundary, and 
mo is the magnetization value at which the droplet shape changes from a rounded 
square to a ring taking advantage of the periodic boundary conditions. 

The obvious qualitative differences in the shapes of i?gp|j(m) and Rili^g{fn), see 

Fig. 3, are the slopes at m = ±m* and the horizontal segment in Ril^j^gim) for 

\m\ < mo. Apparently, the absense of a threshold value tuq in i?gp|^(m) is due to a 
special structure of typical configuration in the spherical model with a fixed value 
of the magnetization. These configurations always take advantage of the periodic 
boundary conditions, and hence they are of the ring shape (as opposed to a round 
droplet shape). The absence of the horizontal segment i?gpj^(m) is due to the diffuse 
nature of the boundary between the "+" and "— " phases. Increasing the amount 
of "+" or "— " phase inside the ring makes the interphase boundary steeper, and, 
hence, makes the set of corresponding configurations less probable. 

Another well-known model, which is often attributed to the mean-field class, 
is the Ising model on a Cayley tree. The explicit expression for the rate function 
i?Cayiey(a;) of the maguetizatiou is not known, nevertheless, it is possible to compute 
the values of this function numerically. It was shown in the paper ^ that Rcayleyi-^J 
is always convex and vanishes at exactly one point where it has strictly positive 
second derivative. Without any doubt the same conclusion is valid for many other 
models on Cayley trees. If so, then one should conclude that ferromagnets on trees 
do not exhibit genuine critical behaviour, where single zero of the rate function either 
stretches to an interval, or splits into several zeroes for sufficiently low temperatures. 
Rather, models on trees are some kind of bundles of (non-critical) one- dimensional 
chains. Distributions of order parameters in tree models always concentrate at a 
single point (the rate-function zero) even for low temperatures and in the absence 
of symmetry breaking perturbations. Non-zero values of order parameters appear 
only as a result of explicit breaking of symmetry by a field applied to an abnormally 
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large number of boundary sites. 

Thus, thermodynamic systems can be classified according to the qualitative be- 
haviour of their large-deviation probabilities. One cap put forward a hypothesis 
that several large classes can be described as follows. Rate functions, R{x), of 
macroscopic observables (order parameters) in non-critical systems are always con- 
vex and vanish (reach zero value) at exactly one point. For sufficiently low tem- 
peratures (strong correlations) the rate functions of (long-range) mean-field models 
are not convex and vanish at a finite number of points. Rate functions of dis- 
crete finite- dimensional models exhibiting critical behaviour are convex and (when 
the temperature is low enough) they vanish on a certain interval [— m*,m*]. If 
[a,b] n [-m*,m*] ^ 0, then 

InPr [mjv G [a,b]] n"^'^ min R'^^^x). 

x£[a,f)] 

Rate functions R{x) of continuous finite- dimensional models are similar to those of 
discrete models. However, if [a,b] fl [—m*,m*] ^ 0, then 

lnPr[mAf G [a, 6]] ~ -nf'^'^ min R^'^^{x). 

a;e[a,6] 

The properties of typical configurations realising large-deviations are even more 
intriguing than the probabilities of these events. Hopefully those will be outlined in 
a consecutive publication. 

Thus, large-deviation probabilities in short-range finite-dimensional systems are 
sensitive to the dimensionality of the lattice. In fact there also exists a dimension 
dependence of another kind. If we look at the magnetization of a subdomain D^^ 
of the entire system then, as it turns out, the behaviour of the corresponding 
large-deviation probabilities is very sensitive to the shape or/and dimensionality of 
Dfc. It is this feature that is the main focuse of the present work. 

The rest of the paper is organized as follows. Section 2 contains the definition 
of the spherical model and statements of the main results. Large deviations of 
the magnetization of various domains arc studied in Section 3. Subsections 3.1 
and 3.2 are devoted to the large-deviation probabilities for the magnetization of 
[d — l)-dimensional layers and [d — 2)-dimensional rods, respectively. Subsection 
3.3 contains derivation of analogous properties for the magnetization of Kadanoff 
blocks. The results of the paper are discussed in Section 4. 

2 Definition of the Model and Main Results. 

Let Ta^ be a rf-dimensional square lattice with nodes j = (ji, ■ ■ ■ , Jd), where jv G Z 
for 1/ = 1, 2, . . . , (i. Consider the sequence of cubes 

^^n = {j e Z'^ : 1 < J, < n, 1/ = 1, 2, . . . , 0?} 

where a random variable G (spin) is attached to each node j G The 
mutual dependence (interaction) of these spins is described by the Hamiltonian 

//„(a) = -J^a,a,-, J>0; (2) 
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where the summation runs over all pairs of nearest neighbours 

d 



u=l 

We assume that the periodic boundary conditions are imposed (for instance, the 
nodes (1,^25 • • • ,jd) and {n,j2, . . . ,jd) are also nearest neighbours). 

The joint probability distribution of the random variables {cxj : j E Qn} is 
specified by the density 

P iWj ■■ J e a,,}) = ^ exp [-(3HrXa)] , (3) 

with respect to the "spherical" a priori measure 

/^nW = 5[ ^ a|-iV| n daj, (4) 

where G„ is the partition function (normalizing factor), = n'^ = \Qn\, is 
the Dirac delta function, and rijen,! d-aj is the Lebesgue measure on (R^,i3(R^)). 
Equations ([2]), ([3]), and (jl]) define the spherical model of a ferromagnet [Tj. 
Consider the sequence of Kadanoff blocks 

Bk = {j enr,:l<j^<k, z/ = 1, 2, . . . , d}, 

where k = k{n) is a non- decreasing function of n (possibly a constant). The nor- 
malized total spin (magnetization) of a block Bk is given by 

Sb = \Bk\ ^ ^ o-j. 
In the present paper we study the distribution densities 

^ +00 

(t)n,B{rn) = — j . . . j 5{Sb - m) exp [-/3if„(o-)] /i„(rfo-) (5) 

" — oo 

of the random variables Sb- In particular, we would like to investigate the asymp- 
totic behaviour of (t)n,B{fn) as n and \Bk\ tend to infinity. 

To carry out the integration in Eq. ([5]) we need to know the spectral properties 
of the symmetric matrix Cn associated with the Hamiltonian ([2]) 

Hn{a) = -J C\fa,a,. 

For the periodic boundary conditions the eigenvalues and eigenvectors of Cn are 
given by 
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where 



1 



Vi f 



27ra-l)(4-l) , . 27ra-l)(4-l) 
cos h sm 



n 



n 



(7) 



The components can be written down in the following compact form 

(2 



Vi 



■ COS 



27r(j.-l)(4-l) TT 

4 



n 



which we will use throughout the paper, although in all calculations we actually 
worked with Eq. (^^. 

No paper on the spherical model on a square lattice Z'^ can be written without 
a reference to the Watson function 



WJz) 



(2vr) 



dui . . . dujd 
-2 - E^!=i cos c;,. 



The integral of this function also appears quite frequently in various expressions and 
certainly deserves a special symbol 



LAz) = 



{2'kY J-n 



dui . . . duJd In 2; — ^ cos uj^ 



(9) 



The main results of the present paper can be stated as follows. Consider the 
Kadanoff blocks Bk and also the following subdomains of the cubes 

Lfc = {j G a„ : 1 < ji < A;}, 

i?fc = {j G fi„ : 1 < ]^ <k, u = l, 2}, 

which we call layers and rods. Let 1 ^ A; ^ n, then the large-n asymptotics of the 
large-deviation probabilities for the magnetization of layers, rods, and blocks are 
given by 

logPr[mL G [a,b]] ~ —n'^~'^ min R^l\x), 



log Pr [niR G [a, h]] 



n . ^(2) 



logn x&[a,b] 



min [x) 



logPr [niB e [a, 6]] ~ -/c"* ^ min R%\x), 

where explicit expressions for the rate functions are given by Eqs. (1231) . (I3T!) . and 
(l36l) . respectively. 



3 Magnetization of Kadanoff blocks. 

To find the distribution density of the magnetization of a Kadanoff block 5, see Eq. 
([5]), we have to calculate the following integral 
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The integration over aj, j G f2n is carried out using the standard technique after 
Berhn and Kac, see Namely, first a new set of integration variables / G f2„ is 
introduced via = J^if^Un ^j^iVh where {Vj-,;}jgn„ = V^, / G fi„ are the eigenvectors 
of the matrix Cn associated with the Hamiltonian Hn{a), see Eqs. ([6]) and ([7]). Next, 
the delta functions are replaced by their integral representations 



1 f+co 

6{a) = — / e^™c?r. 



the integration order is exchanged, and the integration over the variables yi, I & Qn 
is performed. Note that the integration order can be switched only if the quadratic 
form in the argument of the exponential function is negatively defined. This can be 
achieved by a shift of the integration contour for r. One obtains then 



e 



1 



ICO + TQ OO 



dr e' 



dx e 



TT 



-ZOO + TO 



exp 



4(r - (3JXj) 



where 7j = \B\~'^Y,ieB'^j,U cind tq is the shift of integration contour mentioned 
above. For the periodic boundary conditions a straightforward calculation yields 
the following expression for the coefficients 7^, j = (ji, 72, ■ ■ ■ ,jd)- 



1^1 (2/^) n sin[.0.-l)/n] 



cos 



7r(fc-l)(j.-l) TT 
4 



n 



:i2) 



On integrating over the variable x in Eq. ( ITTI) and on introducing a new integra- 
tion variable z via t = /3Jz one arrives at 



On. 



B 



27ri 




dz 



-ioo+zo vSW(^-A(l,...,l), 

where the following notations have been introduced 



■.exp[N(3Mniz,m)], (13) 



E{z) 



E 



7 



z — Xi 



(14) 



2f3JN 



(15) 



ien„\(i,...,i) 

The large-n asymptotic expansion for the remaining integral over z can be derived 
using the saddle-point method. This is quite straightforward when /3 < /3c, where 
/3c = Wd{d)/{2J) is the inverse critical temperature of the spherical model, see [1], 
and Wd{d) is the Watson function at z = d. However, if 



o[n 



some extra 

efforts are required in the low-temperature region /3 > /3c. Difficulties arise because 
the sequence of saddle points {2*} (that is, the sequence of minimum points of the 
function $„(2;,m) on {d, 00)) converges to z = d, as n ^ 00. That is, the sequence 
of saddle points z* approaches the maximal singularity, C,n, of $„(2;, m): z*—C^ni 0, 
as n — )■ 00. Therefore, the standard version of the saddle-point method can not be 
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applied, and the integral ( IT3|) needs a special investigation in the low-temperature 
region. What we actually have to do is to introduce a new integration variable ( via 
z = d + n~^(, where the exponent p must be chosen in such a way that the sequence 
of function $n(C) = ^n{d + n~P(,m) has a "conventional" saddle-point landscape 
in the limit n ^ oo. 

The main goal of the present paper is investigation of large deviation probabilities 
for the magnetization of cubic Kadanoff blocks B. However, it turns out that the 
large-deviation probabilities for the magnetization of layers 

Lk = {j enn-l< Jl < k} 

and rods 

Rk = {j enn:l<Ju<k, z/ = l,2} 

differ qualitatively from those for magnetization of cubic blocks B. Therefore, we 
begin from the simplest (as it happens to be) case of layers L, after that we consider 
large deviations for rods R, and finally we turn our attention to the case of Kadanoff 
blocks B. 

3.1 Large deviations of the magnetization of layers L. 

If the boundary conditions are periodic in all dimensions and 

Lk = {j eQn:l< Jl < k}, (16) 

then a straightforward calculation yields the following expression for the coefficients 
7j = \Lk\-''j:ieL,V,,i (see Eq. (HI])): 



A/2n2 ^ sinWrkiji — l)/n] 

^3 = ~'Tr\ ^TT- TT7T 

Lfc sm7r(ji-l)/n 



7r(A;-l)(ji-l) -k' 
n 4 



J u=2 



Hence, in the case of layers the multiple sum in the expression for '^{z) (see 
Eq. ( IT^ ) reduces to a single one and can be calculated exactly using the method 
described in One obtains 



\Lk\T.{z) = — 



d 



^ 2 {l-x1){x^-xl) 

k{x2 — Xi) (X2 — 1) 



where Xi^2 = I + z — d ^ y {z — d){2 + z — d) . 

Now we are going to locate the maximal singularity of the function $„(-2;, m) given 
by Eq. (flSjl . That will give us a hint how one should rescale the integration variable 
z in order to preserve the saddle-point profile of $„(z,m) in the limit n — )■ oo. Let 
k = an? , where 7 G (0; 1) and a > 0. It is obvious from Eq. that singularities 
of the function Tj{z) — simple poles — are at the points 

z G : 7(j,i,...,i) 7^ 0' J = 1, • • • ,^}- 

To locate zeroes of this function we rescale the variable z according to z = d + (n~'^, 
where C is a new independent variable. The rescaling yields as n — )■ 00: 
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\Lk\T.{d + Cn~P) = 

{ n7C + 0(n3/2p/fc), 

l_exp(-«y2C) 



^27 



n 

T 



1 



, A;n//i(C) + 0{k^ 



for < p < 27, C > 
+ 0(1), forp = 27, C>0 

for 27 < p < 2, C > 
for p = 2, 



(19) 



where 



MO 



2Ctanhv'C/2, 



for C > 0, 

-v/^tan for C < 0. 

Hence, the function T,{z) vanishes at the points 



Zf = d- l-K'il + 2fYn-' + 0{kn 



/ = 0,1,2 



(20) 



Obviously, the integrand in Eq. f|T3|) is an analytic function for z > A(i^...^i). 
The function Ti{z) has a singularity at the point z = A(i_,,,^i). However, for the 
periodic boundary conditions the largest eigenvalue of the interaction matrix Cn is 
non-degenerate, and the singularity of S(z) is cancelled by the multiplier 2 — A(i ...^1). 
Thus, the point z = A(i^... 1) is a removable singularity of the integrand. According 
to Eq. fl20l) the function Il(z) has a simple zero on the interval (A(2,i,...,i); A(i^..._i)). 
Consequently, the maximal singularity of the integrand in Eq. ( fT3l) is at the point 
z = Sn = d — ^TT^n^'^ + 0{kn~^). The obtained asymptotic expansion for the 
location of the maximal singularity suggests that the change of integration variable 
z = d + ^n'"^ might preserve the saddle-point profile of the function $„(z,m) in the 
limit n — )• 00. 

Note now that for any ( > — ^vr^ we have 



$„((i + (n ^, m) 



, Ld{d) _2 



c 1- 



(3 



m^h{C) 



+ 0{kn 



as n — )■ 00, where Ld{d) is given by Eq. ([9]). Hence, the function $„(-2;, m) attains its 
minimum on the interval (s„; 00) at the point 2* = c/ + n^'^C^*{m) + 0{kn~^), where 
C*{m) is the maximal solution of 



/3 



-|7r2 unless m 



0. The point z = 2;,* is the saddle point of the 



Obviously C*("^) > 
integrand in Eq. ( lT3l) . 

In the scale z = d + (n~'^ the saddle point C*(^) does not approach the maximal 



singularity Sn = d 



1^2 



TT n 



+ 0{kn ) of the function $„((i + (n 



m] 



Hence the 



function $n(C) = ^n{d + C'^ ^) ^) has a conventional saddle-point landscape in the 
limit n — 7- 00. Therefore in the scale z = d + (n'"^ one can find the asymptotic 
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expansion for the integral Eq. ( IT3|) using the standard saddle-point method and 
obtain 



Qn,L = exp -iV/(/3) + n'^-'pj 



where 

= i log(/3 J/vr) - pjd + \Ld{d) (21) 

is the limiting free energy per spin of the spherical model for /3 > /3c- The parti- 
tion function of the spherical model with the periodic boundary conditions in the 
canonical ensemble is given by 



e„ = exp{-iV/(/3) + 0(logn)} 



(22) 



Thus, we obtain the following asymptotics for the distribution density of the mag- 
netization of the layer L, see Eq. 

</)„,i(m) = exp + 0(fcn'^-3)} , 

as n — )■ 00, where the rate function is given by 



Rf{m) = -(3 J 



(23) 



(2) 

The function R\ [m) has a cusp at m = and vanishes at the points m = 
Pel P — the equilibrium values of the magnetization within the spherical 



model with periodic boundary conditions, see Fig. 3. The nature of the cusp at 
m = is obvious. A magnetization value m can be achieved by deforming either 
equilibrium (+)-phase or (— )-phase. For negative (positive) values of m a defor- 
mation of minus (plus) phase is the easiest way to achieve the desired value of 
magnetization. As m the required deformations become more and more costly, 
because we are getting further and further away from the equilibrium value —m*. 
However, if we keep increasing the value of m deformations become less costly as 
soon as we cross the point m = 0, because we begin to approach the equilibrium 
value m*. 



3.2 Large deviations of the magnetization of rods R. 

In this subsection we investigate large deviations of the magnetization of {d — 2)- 
dimensional domains (rods) 

Rk = {j enn:l< 3u <k, z/ = 1, 2} . 

As in the previous subsection the large-deviation probabilities are determined by 
the relative location of the two largest eigenvalues of the interaction matrix and the 
maximal zero of the function Ti{z). However in the case of rods Rk the distances 
between the eigenvalues and zeroes are qualitatively different from what we found 
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Figure 4: The rate functions for the magnetisation of the layers L (left) and the 
Kadanoff blocks B (right) at /3 = \(3c. 



in the case of layers L^. As a result, the probabilities of large deviations are also 
qualitatively different in these two cases. 

In the case of rods the coefficients 7^- = |-Rfc|^^ HieRk ^'.^ given by 



7i 



„ sin[7rA;(j^ - l)/n] 

I I —r — . . , , COS 



|i?fc|n'^/2-2 ^11^ sin[7r(j^ - l)/n] 



7r(A;-l)(j,-l) 



. K=3 



Hence the multiple sum in the expression for S(-2) (see Eq. ( |T4l) ) reduces to the 
double one 



S(z) 



d-4 



E 

h,32= 



sin^[7rfc(ji — 1)1 n] sin^[7rA;(j2 — 1)/''^] 
1 sin^[7r(ji - l)/n] sin^[7r(j2 - l)/n] 



z — d + 2 — cos[27r(ji — l)/n] — cos[27r(j2 — 1)/''^] 

The function S(z) has simple poles at the points z = d and z = d — 1 + cos(27r/n) 
and a simple zero s* somewhere between the poles. Our objective now is to find an 
expression for ^(2;) convenient enough for locating s*. 

In the case of {d — l)-dimensional layers the expression for Tj{z) contained only 
a single sum, which we managed to calculate exactly. For {d — 2)-dimensional rods 
we do not have a luxury of working with an exact closed-form expression for S(z), 
nevertheless, fortunately, simplifications arise for a different reason. As is often the 
case for double sums of this type S((i + (n'"^) = 0{n'^~'^k^ log n), for > 0. Under 
the same rescaling the term corresponding to ji = j2 = (the contribution from 
the maximal eigenvalue of the interaction matrix) is given by n'^~'^k^/C. Denote 
now Tj'{z) the same sum as in Eq. ( !T4|) but with the contribution from the maximal 
eigenvalue omitted. Below we will argue that the main asymptotics of S'((i+C^~^) is 
given by n'^~'^k'^\ogn(j{Q, where cr(C) is a positive analytic function for > —2tt'^. 
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Therefore, in order to find s*, we have to solve the equation a{Qlogn — — ^ ^. 
Since the function a{Q is regular at ^ = 0, the solution is given by 



1 



and hence 



(7(0) logn 
1 



+ 



^log^ n J ' 



as n — >■ 00. 



(T(0)n2 logn' 

Let k = arC , where a > and 7 e (0; 1). For z = d + C,''^~^i where C > ^^'^ 
< p < 2, one has as n ^ 00 



exp 



2Cn 



l-p/2 



where 



2-K 



For z — d-\- (n~'^ one obtains 



1 — cos(A;a;i,) 



cos UJ; 



X + m — Yl cos Uk 

K=l 



(25) 



n 



,d-2 



where 



and 



HO - E 



{/2(Cn-2) + k' [6,{Q + ^2(0 + 0(n-(^-^)/2)] } , 
2 1 



J=-0O 



y2(C + 2^ exp y2(C + 2^ - 1 ' 



1 



TT >/0 



duj 



(26) 



Note that A(i^,,,^i) — A(2,i,...,i) ~ 27r^n^^, hence, to locate s* we have to consider the 
scale z = d + C^~^ or even a finer one. In any rougher scale the poles and zeroes of 
^{z) merge as n — > 00, and the function ^n{d + C^"'') "^) does not have a required 
saddle-point landscape. 

Isolating the singularity of the integral l2{C,n~'^) one can write it down in the 
following form 



k' 



dujiduj2 



Cn-2a-7) + i(a;2 + a;i) 



fc'r«(C), 



where the function r^n\C) (the regular part of the integral) is analytic in the unit 
circle {C, E C : \C, \ < 1}. Asn— )-oo the sequence r^\C,) converges to some limiting 
function uniformly over the unit circle. Calculating the double integral one obtains 



2tt 



log^^ + log(ic2 + Cn-^(^--)) 



+ kV:\C). (27) 
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Separating the singularities of 5i(C) at ^ = one obtains 

5i(C) = ^-^ + r(^nC), (28) 

where r^'^^Q is analytic for Re C > — 27r^. Analogous separation of singularities at 
C = yields 



^2(0 = - / 
71 J 



duj 



du 



=5 + '■3(C) 



^ + il^ + ?,(0. (29) 



where the functions r^lQ and r^lQ are analytic in the unit circle. Summarizing 
Eqs. ([27D, (I2S]), and di]) one arrives at 

■1-7 



logn + C ^ + r„(C) 



(30) 



where r„(C) are analytic and bounded uniformly over n in the unit circle. Hence, 
the maximal zero of the function Ti{z) is at the point 

1 



vr 

^n = d- ip. 

1 — 7 log n 



+ 



[nlogn) 



The expression for s* suggests that the change of variable z = d + ({n'^logn)^^ 
is likely to convert Eq. f|T3|) into an integral convenient for application of the saddle- 
point method. Substituting z = d + ({n"^ logra)"^ in the expression for $„(2;) we see 
that this is indeed the case. One obtains 

c \ 



d + 



71^ log n 



d- 



m 

2/3J 



+ 



n''-'l3J 
\ogn 



Tim 



1 - 7 + ttC 



-1 



O 



n 



d-2 



log n 



Hence, the sequence of relevant saddle points is given by 



ra^ log n 



where 



TT 



+ 



1 



n'^ log n , 



VVl-/3c//3 

Evaluating the integral in Eq. (IT5]1 using the saddle-point method one obtains 



<$)n,r = exp 



NfW) 



m\ 



log n 1 — 7 

Taking into account Eq. fl22]) one arrives at 



n 



^ log^ n ^ 



where 



\R[m) = exp 



i?g)(m) 



iyjd—2 



n 



d~2 



logn 



log n 



n(3J 



m\ 



(31) 
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3.3 Large deviations in KadanofF blocks. 

In this section we investigate the probabihties of large deviations for the total spin 
of the Kadanoff blocks 

Sfc = {j e fi„ : 1 < < A;, u = l,...,d}. 

Qualitative behaviour of the corresponding large-deviation probabilities is the same 
for any d > 3. Therefore, to avoid unnecessary technical complications we consider 
only the case d = 3. The coefficients 7^ = |-Bfc|~^ ^leBt. Vj,i ^^^^ given by 



_ 2^/^ sin[7rfc(j. - l)/n] 

\B.\Ji, sin[vr(j. - l)/n] 



■7r(A;-l)(j.-l) TT 



n 



(32) 



and 



1 



j^(^) = ^-3 A sin^[7rA:(j.-l)/n] 

^ sin2[7r(j,-l)/n] E cos[27r(j« - l)/n] 



Jl J2,J3 = 



K=l,2,3 



The function T,{z) has simple poles at the points z = d and z = d — 1 + cos(27r/n) 
and a simple zero s* in the interval {d—l + cos(27r /n);d). Our goal now is to obtain 
an expression for Ti{z) convenient enough for locating s*. 

Let k = arC , with a > and 7 G (0; 1). For z = d + (n'"^, C > 0, one has as 
n — )■ 00 



E(d + Cn-') = UCn-') + - 6,{C) + HO + h{C) + 0{n 



-(l-7)/2^ 



n 



where Isi^x) is given by Eq. (125|) . 



1 



U=-oo J2C + 47r2(j2 + /2) exp + 47r2(j2 + /2) _ 1 



52(c) = ^ E / 



j=-oo 



and 



^3(0 = ^ duidu: 



V2C + 4:7r'^f + exp v/2CT4pj^Tu^ - 1 ' 
2 1 



'2C + + exp y^2C + + - 1 
The large-n asymptotics of Is{(n~P) for p > 27 is given by 

where 



2 ^1 
^3 = ^ // / duJiduJiduj'i W — —3 



1 



(33) 



(34) 



(35) 
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Separating the singularities of 6u{(), = 1, 2, 3 at C = (cf., Eqs. (!27|) -( l29|) ) one 
obtains 

E '^.(C) = r^ + r(c), 

u=l,2,3 

where the function r{() is analytic in the unit ball {( E C : \(\ < 1}. Hence, the 
maximal zero of the function T,{z) is at the point 



Z = Sn ^ d 



k 



The location of the maximal zero gives us a hint that the asymptotic expansion of 
the integral ( [T3|) can be found using the saddle-point method after a prior rescaling 
of the integration variable z via z = d + (n~^'^'^ . The rescaling yields 



L{d) 
2PJ 



+ kn 



P ^3 + 



+ o{kn~^). 



Hence, for (3 > f3c the sequence of relevant saddle points is given by 

z*^ = d + Ckn'^ + o{kn~^), 

where 



W3 



1 - /3c//3 



Evaluation of the integral f|T3|) using the saddle-point method yields 



<S)n,b = exp 



NfiP) - k^ (\m\ - Vl-/3e//3)' + oik) 



Taking into account Eq. (!22|) . one obtains the following asymptotic formula for the 
probability density of magnetization 



,B = exp 



where 



W3 



-kR's\m) + o{k)] , 
[\m\-,Jl-/3J/3y. 



(36) 



4 Discussion and concluding remarks. 

In the present paper we have investigated large-deviation probabilities for the mag- 
netization of various domains -Dfc(c Qn) within the spherical model of a ferromagnet 
defined in cubes f2„. We have shown that these probabilities are very sensitive to 
the shape and/or dimensionality of the domains D^.. 

An appealing feature of the large-deviation theory is that often the behavior of 
large-deviation probabilities admits a simple intuitively clear explanation, see, e.g., 
|8] . This is also the case for the results obtained in the present paper. In Section 
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3.1 we have shown that in the case of layers Lk the large-deviation probabilities 
decay exponentially with n'^~^, which is similar to the asymptotics of the analogous 
probabilities for the magnetization of the entire cube see [B]. The reason for such 
similarity becomes clear when we investigate the properties of typical configurations 
realizing these deviations. It turns out that the easiest way to obtain a desired value 
of magnetization in the layer Lk is to deform the configuration of random variables 
in the entire cube Qn- That is, a large-deviation of the magnetization in a layer 
Lk leads to a large-deviation of the magnetization in the entire cube Qn- To put it 
another way, a substantial deformation of the configuration in a layer L^ spreads 
over the entire cube fin- 

In the case of rods Rk the large-deviation probabilities are modified by log n cor- 
rections, namely, they decay exponentially with n'^~^/logn. Investigation of typical 
configurations realizing these large deviations shows that a substantial deformation 
of the configuration in a rod Rk does not spread over the entire cube but it 
spreads over a domain with the linear size of the order n/ logra. Note that if k satis- 
fies the bound y/nk' <^ <^ n, then the number of random variables in a rod Rk is 
greater than that in a layer Lk'. Nevertheless, large deviations of the magnetization 
in the rod are more likely than large deviations in the layer simply because of a 
more compact arrangement of the random variables in the rod. 

In the case of Kadanoff blocks Bk the large-deviation probabilities decay expo- 
nentially with k'^~'^. Even a substantial deformation of the configuration in a block 
Bk does not spread very far from the block. Essentially it remains localized in a 
domain of the linear size 0{k), that is, of the same extent as the size of the Kadanoff 
block itself. The number of random variables in a block Bk can be greater than the 
number of random variables in a rod Rk' (or a layer Lk"), nevertheless, a large de- 
viation of the magnetization of the block is much more likely than a large deviation 
in the rod or the layer. Accordingly, a uniform deformation of the configuration in 
a layer Lk" or a rod Rk' spreads much further over the cube fin than a deformation 
of a Kadanoff block Bk- 

In order to find the properties of typical configurations realizing large deviations 
in layers, rods, and blocks one has to continue the calculation of the present paper 
one step further and derive the conditional distributions of random variables aj 
given a desired large deviation. These calculations will be published elsewhere. 
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